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We investigate numerically dynamical aspects of the discrete nonlinear Schrodinger equation (DNLS) in finite 
lattices. We start from a periodic chain with nearest neighbor interaction and insert distant links that connect 
random pairs of sites across the lattice. For a localized initial condition we focus on the time averaged probability 
at the initial site analyzed as a function of the degree of complexity of the lattice. We find that selftrapping occurs 
at increasingly larger values of the nonlinearity parameter as the lattice connectivity increases while close to the 
fully connected network localization becomes more preferred. The nonlinear network dynamical behavior is 
understood as a competition of two lattice limits, viz. that of pure nearest neighbor interaction as well as the 
mean field limit where all sites are mutually coupled. 



I. INTRODUCTION 

The aim of this work is to investigate what happens when 
two distinct aspects of complex systems coexist simultane- 
ously. The first feature is that of selftrapping, viz. a generic 
dynamical property of nonlinear lattices that leads to local- 
ization in a translationally invariant system. The second is 
the geometrical feature that emerges from a small world net- 
work (SWN). In order to probe the former we use the Dis- 
crete nonlinear Schrodinger (DNLS) equation that is an ubiq- 
uitous equation arising in a variety of contexts in condensed 
matter [d, HI HI, m optics 01, in Bose-Einstein condensation 
JH], etc. For the second feature we use Newman- Watts Small- 
world networks (NWSWN's) where random links are inserted 
linkingpairs of the original nearest neighbor coupled lattice 
sites ESI- 

The DNLS equation is given by: 



i-^- = V(C„_i + C n+ i)-%\C„\ 2 C n 



(1) 



where the parameters % and V respectively introduce anhar- 
monicity and dispersion; hereafter we set V = 1. The C n (t) 
are complex amplitudes with C„| 2 being the probability for 
finding the excitation at site n at time t. We take the total 
probability normalized to unity £ C„| 2 = 1. We may express 
the DNLS equation also in matrix form yfl: 



/C=MC-x£>(|C| 2 )C 



(2) 



where C = col(C\,C2, • • • ,Cjv) is an N-component col- 
umn vector, each component representing the complex ro- 
tating wave amplitude of a particular oscillator, D(|C| 2 ) = 
diag(\Ci,C2, ■ ■ ■ ,Cff\) is a diagonal N x N matrix represent- 
ing the nonlinear term which introduces anharmonicity and 
M = [M mn ] is a symmetric N x N matrix, the elements of 
which represent linear interactions between oscillators. 

One of the most interesting features of DNLS is that it leads 
to self-trapping, viz. while for small values of % the probabil- 
ity wave is delocalized, this does not hold for larger nonlinear- 
ity values yfl.The self-trapping properties of DNLS for finite 
chains of sites have been studied extensively in Ref @]. But 
single lattice models are often insufficient to describe more 



complex systems encountered in nature. In general such sys- 
tems are characterized by graphs with bonds connecting sites 
with a wide distribution of mutual distances. Examples can be 
found in various fields, ranging from physics or biology to so- 
cial science and computers. This structural complexity gives 
rise to interesting statistical properties. These systems are 
known as small-world networks, i.e. networks whose bonds 
posses a high degree of local clustering, and at the same time 
a short average path length connecting the nodes of the system 

HEzLi. 

In the present paper we extend the work of Ref. U to struc- 
turally complex networks. Specifically, we focus on the de- 
pendence of localization properties of DNLS on the degree of 
complexity of the network; the problem is set as follows: We 
begin with N lattice sites that form a one dimensional ring. 
For relatively large values of N the DNLS dynamics on this 
ring is not different form the one in very long one dimen- 
sional lattices @]. Subsequently, we destroy the periodicity 
of the lattice by randomly inserting connections between dis- 
tant sites, as seen ion Ref fioll . As an initial condition we use 
that which places the probability initially on one site of the 
system. The main quantity studied throughout this paper is 
the long time-averaged probability Po(t), viz. the occupation 
probability of the initially populated site: 



(3) 



Physically, we are placing a quantum particle at a given 
lattice site and follow its time evolution in the semiclassical 
limit. This study is related to work done on DNLS in disor- 
dered lattices 



t udy is relatei 



II. SELFTRAPPING IN A COMPLEX LATTICE 

Before we address the dynamics of the DNLS equation in 
a small world network we need first to establish two specific 
limits. In order to form a small world network, we start from 
a periodic ring of sites and add randomly (using a uniform 
distribution) B bonds linking distant sites. For B = only 
first neighbors are connected; this is the nearest neighbor limit 
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(NN). On the other hand, the system is fully connected for the 
maximum number of possible distant connections, viz. for 
Bmax = N x (N — 1 ) /2 — N; we refer to this as the mean field 
(MF) limit. In the NN limit we have the DNLS equation on a 
regular lattice with periodic boundary conditions. In this case, 
when nonlinearity is zero, an initially localized state spreads 
throughout the lattice and after long times the probability to 
find it at a given site will be of the order of 1 /N. Increasing of 
the nonlinearity parameter % leads to breather formation sig- 
naled by selftrapping for values in the ra nge % w 3 — 5 121 . 
In the MF limit, on the other hand, the network is fully con- 
nected and this has ramifications for the DNLS selftrapping 
properties. We note that the pure linear chain on a small world 
network has been investigated in Ref ifioll . It was found that 
for B small partial localization of the probability at the initial 
site occurs. 

In order to study the general case of the DNLS on a SWN 
we use the following numerical procedures: We start from 
the chain with nearest neighbor couplings and add randomly 
bonds to distant sites. We place initially all probability at a 
given site (referred to as the "zeroth" site) and for each lat- 
tice we follow the evolution of the probability amplitudes as 
a function of time. We evaluate Po(t), viz. the probability at 
the initial site, for long times (T max = 10 5 ) and time-average 
this this probability in order to obtain < Pq >. This numer- 
ical experiment is repeated for R distinct lattice realizations 
and, finally, we obtain the averaged self-probability P over 
these realizations. We present results for N = 100 while the 
ensemble averaging is done over R = 500 realizations. Re- 
peating the above procedures for different values of the non- 
linearity parameter % results in the behavior seen in Fig. (Q]). 
We observe that the well known selftrapping transition of the 
finite, extended lattice (B = 0) occurring for x — 3.5 H] shifts 
to larger values of the nonlinearity parameter as we increase 
the distant bond number B. In other words, the onset for the 
breather regime close to the NN limit becomes delayed by 
the presence of a small number of distant bonds. This be- 
havior is understood since, in this limit, the addition of few 
distant bonds, even though adds disorder to the system, it nev- 
ertheless, increases the connectivity of the lattice and, as a 
result, the amount of nonlinearity necessary for selftrapping 
is higher. 



III. FROM THE NEAREST NEIGHBOR TO THE MEAN 
FIELD LIMIT 

In order to explore the whole range of complex lattices from 
B = to B = B max we need to focus first on the dynamics 
at the fully coupled MF limit where each site is connect to 
every other site. In the purely linear case (% = 0) the nonlinear 
eigenvalue problem of Eq. (f2j) reduces to 



MC = XC det[M -XI] =0 



(4) 
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where C is the n-component column vector and V is the inter- 
action matrix. This system is solved easily leading to N — 1 
degenerate localized eigenvalues with Xi oc = — 1 and a single, 
delocalized (Perron-Frobenius) eigenvalue for X^el = N — 1 . 



FIG. 1 : Time and realization averaged probability Pas a function of 
the nonlinearity parameter % for different dilute SWN lattices. Self- 
trapping occurs at increasingly larger values of X as a function of 
B. 



In the large-Af limit thus, initial conditions remain fully local- 
ized except for the special case where the particle is placed 
equally on all sites; in the latter case it remains delocalized. 
We thus find that in the linear case while the NN limit favors 
derealization the MF limit promotes localization! This fea- 
ture will thus dominate the behavior of DNLS in the large-Z? 
limit. 

In order to study the dynamics of DNLS towards the MF 
limit we employ two distinct approaches. In the first we sim- 
ply change the number of bonds B to B max while in second 
we add all B max bonds in the system from the beginning while 
assigning to each one a given strength Wb- Subsequently we 
change Wb from zero (no small world links) to one (MF limit). 



A. Variable distant bond number 

We construct our lattice by starting from the NN limit 
where B = and add randomly bonds until we reach the MF 
limit where B = B max (Fig. (2)). We use N = 100 and thus 
Bmax = 4850. The linear case (% = 0) for relatively small B 
has been addressed by Blumen et al. [flOfl where is noted that 
the addition of a few distant bonds, while primarily delocal- 
izes the probability, it nevertheless causes small, partial local- 
ization, a phenomenon that could be connected with Anderson 
localization fl4ll . We also observe this behavior, viz. that ap- 
proximately 5% of the total probability remains at the initial 
site (fig. ID in the small-B limit, while the majority of the prob- 
ability wave spreads in the network. However, as the number 
of bonds increases and upon approaching the MF limit a rela- 
tively sharp transition occurs for values of B very close to the 
MF, as can be seen in Fig. ([3^). The averaged probability, 
while being essentially delocalized for all intermediate B val- 
ues, it now collapses back to the initially populated site. As 
a result, when the linear network is close to the MF limit it 
tends to trap the probability in the initally populated site. 
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FIG. 2: Starting from the NN limit where B - 
bonds until we reach the MF limit where B - 



and we add randomly 

- B/nax- 



We now turn the nonlinear term on and observe in the 
regime from B = to B close to B max two competing tenden- 
cies, viz. derealization due to bond addition as well as local- 
ization due to nonlinearity. For small number of bonds nonlin- 
earity induces localization that nevertheless occurs for larger 
values of % as the bond number increases. In the intermedi- 
ate B regime we observe derealization unless the nonlinearity 
parameter values become very large. In this range nonlinear- 
ity competes strongly with the derealization induced from 
the distant bonds and a minimum is reached where the sys- 
tem acquires its most delocalized form. Finally, in the large 
bond number regime the MF regime augments the localization 
properties of the lattice leading to selftrapping for smaller val- 
ues of the nonlinearity parameter. This rich behavior is shown 
in Fig|3j where in [3J a ) we see the full range P behavior as a 
function of the bond number for different nonlinearity values 
while in [3{b) and 0c) we observe the small (NN) and large 
(MF) behavior respectively. We note specifically the rapid lo- 
calization feature in the MF limit. 
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FIG. 3: Time and realization averaged probability of occupation of 
the initial site as a function of the bond number B for different non- 
linearity values, (a) Full bond number range, (b) detail of the nearest 
neighbor regime and (c) detail of the mean field regime. 



B. Variable distant bond strength 

We employ additionally a second approach whereby we 
start from the full SWN system assigning however a certain 
weight Wb to all the distant bonds ranging from Wb = (NN 
regime) to Wb = 1 (MF limit) as seen in Fig.©. We find qual- 
itatively similar results as in the previous case as shown in 
Fig-©; the selftrapping transition occurs at distinct values of 



Wb depending on the nonlinearity parameter values. As the 
value of nonlinearity increases the onset of localization occurs 
at smaller values of Wb as expected. If we follow a curve for 
specific % we note that first we enter in the selftrapping regime 
while for larger values of Wb we get an additional "transition" 
stemming from the MF attractor. This further localization oc- 
curs for smaller Wb values as the nonlinearity parameter in- 
creases. 




FIG. 4: Starting from the fully connected SWN system assigning 
however a certain weight W;, to all the distant bonds ranging from 
W b = (NN regime) to W b = 1 (MF limit) with V = 1 . 

We notice that as we gradually increase the W/,'s that the 
probability of the initially pertubated site mark a rapid change 
for values very close to one. 
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FIG. 5: Time and realization averaged probability of occupation of 
the initially populated site as a function of the bond strength Wj, for 
different nonlinearity parameter values. 



IV. CONCLUSIONS 

We investigated dynamical properties of the DNLS equa- 
tion related to selftrapping and the onset of nonlinear localiza- 
tion in complex lattices in the form of small world networks. 
The nonlinearity parameter % and the number of small world 
connections B (or their strength Wb) are the two variables that 
induce localization and derealization respectively. 

We saw that in the fully linear regime (% = 0) the increase 
of the small world link number leads essentially to derealiza- 
tion due to the effective increase of the lattice dimensionality. 
This trend is true except very close to the mean field limit; in 
the latter, exact solution shows that the lattice has complete lo- 
calization tendency (except in the case of a totally delocalized 
initial condition). We note, however, that in the very small B 
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regime, while the system is delocalized, some partial local- 
ization is observed flOll . When the nonlinearity parameter is 
non zero we observe a competition between localization in- 
duced by the latter and derealization induced by the small 
world connections. This competition makes selftrapping oc- 
cur at larger %-values compared to the translationally invariant 
lattice case while for very large nonlinearities the probability 



becomes to a large extent selftrapped. An interesting feature 
appears from the linear MF behavior, viz. close to the large 
B regime the nonlinear localization tendency is augmented by 
the linear MF localization behavior. These features are also 
seen in the fully linked lattice with variable bond strength. 
The DNLS properties on the SWN may have some applica- 
tions in energy propagation and storage in complex networks. 
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